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I. INTRODUCTION 



Unlike ordinary three-dimensional systems, quantum two-dimensional systems of indis- 
tinguishable particles allow for generalized braiding statistics. A celebrated generalization 
of the usual bosonic and fermionic quantum statistics is provided in two dimensions by 
Abelian anyons, for which a phase factor multiplying the scalar wavefunction is associated 
to elementary braiding operations [l ^] 



l/j(Zi, .., Zi, . . . , Zj, .., z n ) 



tp(zi, .., Zj, . . . , Zi, .., Z n ) ■ 



(1) 



Anyons, first studied in 



4j-|6|, were later associated to the physics of the fractional quantum 



Hall effect Abelian anyon statistics of the simplest QH states, at filling factors v = 
1 / (2p + 1) were derived from a microscopic theory |z|: since then, the study of the properties 
of Abelian anyons and the applications to the QHE have been in the following decades subject 
of an intense and continuing interest 12 1 . 

A further generalization of the bosonic and fermionic statistics is represented by 
non-Abelian anyons, described by a multi-component wavefunction ip a (zi, . . . , z n ) (a = 
1,2, . . . , g) which undergoes a linear unitary transformation under the effect of braiding 
cTj which exchanges the particles at the positions Zi and Zi + \ 



(2) 



where p(&i) are g x g dimensional unitary matrices which do not commute among themselves, 

[p((?i)U[p(Vj)}bc ^ [p(Vj)]ab[p(Vi)]bc Q • 

Abelian and non-Abelian anyons respectively correspond to one- dimensional and higher- 
dimensional representations of the braid group: with respect to parastatistics, non-Abelian 
anyons represent the counterpart of the generalization represented by Abelian anyons with 
respect to ordinary Bose and Fermi statistics. Non-Abelian anyons naturally appear in the 



description of a variety of physical phenomena, ranging from the 



scattering of vortices in spontaneously broken gauge theories 



'ractional QHE |3|,ll3[ to the 



14 



gravity 



-1 161], the (2+l)-dimensional 



17H19| and the alternation and interchange of e/4 and e/2 period interference oscil- 



lations in QH heterostructures 



20] 



The investigation of the proprieties of non-Abelian anyons is currently attracting a wide 
attention, both per se for the peculiarities of such strongly correlated states, but also for 
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,heir relevance for the possible implementation of topological quantum computation schemes 
. The reasons for such an interest are related to the fact that, using the braiding rules of 



non-Abelian anyons (such as the Fibonacci anyons [3]), one can implement all the operations 
needed for universal quantum computation 2l| : the efficiency of schemes to perform qubit 
operations has been discussed in 22|, |23|. The main point is that non-Abelian excitations 
are topological, and therefore expected to be rather robust with respect to the decoherence 
caused by the interaction with the environment. Strong evidences of the non-Abelian nature 



of the excitations in the 5/2-state have been found in QH systems [2J, |25 |. 
the recent realization of synthetic gauge fields in systems of ultracold atoms 



urthermore, 



26] opens to 



the possibility of implementing theoretical proposals for the realization and manipulation 



of non-Abelian anyons in such systems 27H30| (see more references in the reviews 31 



32J). Proposals for stabilizing non-Abelian anyons on the edges of fractional quantum Hall 



33 



34[ have been 



states by proximity-coupling to superconductors and/or ferromagnets 
very recently discussed in literature. 

The study of equilibrium properties of two-dimensional systems having non-Abelian ex- 
citations is in general a non trivial and highly interesting task: indeed, the anyonic statistics 
incorporate the effects of interaction in microscopic bosonic or fermionic systems (statistical 
transmutation) so that the determination of thermodynamical properties of non-interacting 
anyons is at least as much as difficult as similar computation in ordinary interacting gas. 
This is a reason for which the investigation of e quil ibrium properties of a free gas of anyons 
called for an huge amount of efforts and work [35[, the other reason of course being that 
the thermodynamics of a system of free anyons is the starting point - paradigmatic for the 
simplicity of the model - for the understanding of the thermodynamics of more complicated 
interacting anyon gas. 

The two-dimensional gas of free Abelian anyons whose wavefunction fulfills hard-core 
wavefunction boundary conditions has been studied by Arovas, Schrieffer, Wilczek, and Zee 
36| in its low-density regime by taking its virial expansion. In particular, they found the 
exact expression for the second virial coefficient, that turns out to be periodic and non- 
analytic as a function of the statistical parameter. Results for higher virial coefficients of 
the free Abelian gas are also available in literature: different approaches have been used, 



including the semiclassical approximation 
references see 



37J and Monte Carlo computations [38] (for more 



39[). Useful results can be found by perturbative expansions in powers of 



i 



the statistical parameter a: exact expressions for the first three terms of the expansions in 



40 



4l| . The second virial 



powers of a are available for each of the first six virial coefficients 
coefficient is the only one presenting - in each of the Bose points - cusps in the statistical 
parameter a i.e. none of the higher virial coefficients have terms at order a 42, 43]. 
Furthermore, a recursive algorithm permits to compute the term in a 2 of all the cluster and 
virial coefficients 142- 



45]. 



The results for the virial coefficients of the free gas of Abelian anyons quoted in the 
previous paragraph are obtained considering a many-body anyonic wavefunction fulfilling 
hard-core boundary conditions, i.e. a wave function which vanishes in correspondence of co- 
incident points in the configuration space of the set of anyons. The generalization obtained 
by removing such an hard-core constraint has been studied for Abelian anyons 46l448| and 
a family of anyon models can be associated to the different boundary conditions of the 
same Hamiltonian. These models are obtained within the frame of the quantum-mechanical 
method of the self-adjoint extensions of the Schrodinger anyonic Hamiltonian. In the fol- 
lowing we will refer to anyons without the constraint of hard-core conditions as "soft-core" 
or "colliding" anyons. The mathematical arguments underlying the possibility of such a 



generalization were discussed in 
anyons was studied in 



47 



48 



46[, and the second virial coefficient of soft-core Abelian 



. The corresponding self-adjoint extensions for the non- Abelian 

45 1| . We stress that it is not easy, in gen- 



anyonic theory have been thoroughly discussed |49 
eral, to extract the parameters of emerging effective (eventually free) anyonic models from 
the microscopic Hamiltonians, and then the introduction of soft-coreness conditions may 
provide useful parameters which have to be fixed via the comparison between the results of 
the anyonic models and the computations done in the underlying microscopic models. 

For non- Abelian anyons, a study of the thermodynamical properties in the lowest Landau 
level of a strong magnetic field has been performed 52|, showing that the virial coefficients 
are independent of the statistics. The theory of non-relativistic matter with non-Abelian 
Chern-Simons gauge interaction in (2 + 1) dimensions was studied adopting a mean field 
approximation in the current- algebra formulation already applied to the Abelian anyons and 



finding a superfluid phase 



53]. 



In comparison with the Abelian case, the thermodynamics of a system of free non-Abelian 
anyons appears to be much harder to study and all the available results are for hard-core 
boundary conditions 54m7| , with - at the best of our knowledge - no results (even for the 
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second virial coefficient) for soft-core non-Abelian anyons. 

The reason of this gap is at least twofold: from one side, for the difficulties, both analytical 
and numerical, in obtaining the finite temperature equation of state for non-Abelian anyons 
(see the discussion in 35]); from another side, because most of the efforts have been focused 
in the last decade on the study of two-dimensional systems which are gapped in the bulk 
and gapless on the edges, as for the states commonly studied for the fractional quantum Hall 
effect, while, on the contrary, the two-dimensional free gas of anyons is gapless. However, 
there is by now a mounting interest in the study of three-dimensional topological insulators, 
systems gapped in the bulk, but having protected conducting gapless states on their edge 



or surface [58j: exotic states can occur at the surface of a three-dimensional topological 
insulator due to an induced energy gap, and a superconducting energy gap leads to a state 
supporting Majorana fermions, providing new possibilities for the realization of topological 
quantum computation. This surging of activity certainly calls for an investigation of the 
finite temperature properties of general gapless topological states on the two-dimensional 
surface of three-dimensional topological insulators and superconductors. 

In this paper we focus on the study of the thermodynamics of an ideal gas of a general class 
of non-Abelian Chern-Simons (NACS) particles in presence of general soft-core boundary 
conditions: explicit results are found for the second virial coefficient. Results for hard-core 
non-Abelian anyons, which is a limiting case of soft-core conditions, are presented too. The 
article is structured as follows: in Section [III we introduce the NACS model studied in the 
paper and, as an introduction to the subsequent discussion, in Section III Al we briefly recall 
the results for an ideal gas of hard-core Abelian anyons and we present in detail the general 
soft-core version of the Abelian anyonic model. The properties of the virial expansion are 
also reviewed and the monotonic behaviour of the second virial coefficient B 2 with respect to 
the statistical parameter is taken in exam as the degree of hard-coreness changes: we observe, 
in particular, that for a narrow range of the soft-core parameter B 2 can be no n- monotonic. 
In Section III Bl we define more precisely the NACS model and we explicitly present the set 
of soft-coreness parameters associated to the most general boundary conditions of the wave- 
functions. In Section II II Al the coefficient B 2 is evaluated for a system of NACS particles 
with hard-core boundary conditions: we compare our results with previous determination of 
this quantity and we make some comments about limit cases. In Section IIII Bl we study B 2 
for non-Abelian anyons when soft-core wavefunction boundary conditions are allowed, with 
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special attention to the case of isotropic boundary conditions. In Section IIVI we summarize 
the virial expansions for the ideal gas of NACS particles endowed with general boundary 
conditions of the wave functions. Our conclusion are discussed in Section [V] Finally the 
Appendices deal with some technical details of the main text and with the energy spectrum 
in the soft-core case. 

II. THE MODEL 

In this Section we introduce the Abelian and non-Abelian models studied in the paper: 
in Section THAI we first briefly remind the well-known results for the thermodynamics of the 
ideal gas of hard-core Abelian anyons. The general soft-core version of the Abelian anyonic 
model is then introduced, and the behaviour of the second virial coefficient is studied as a 
function of the defined hard-coreness parameter. In Section Hi Bl we define the NACS model, 
whose second virial coefficient will be derived and studied in the next Section. 



A. Abelian Anyons 



or a system of identical Abelian anyons has been developed starting 



36[, in which the exact quantum expression for the second virial 



The thermodynamics 
with the seminal paper 
coefficient is derived: 

B h 2 c {2j + 6,T) = -1a| + \5\X 2 T - l -5 2 \ 2 T . (3) 

Eq. 03]) holds for an ideal gas of anyons whose wavefunction fulfills hard-core wavefunction 
boundary conditions. In a = 2j + 5, where a represents the statistical parameter of 
the anyons Q], j is an integer and \5\ < 1. We remind that a =JJ and a = 1 corresponds 
respectively to free two-dimensional spinless bosons and fermions 
thermal wavelength defined as 

2nh 2 N 1/2 



2|. Furthermore At is the 



Ar =fcJ ' (4) 

As discussed in statistical mechanics textbooks, the virial expansion is done in powers of 
pA T (where p is the density and M is the mass of the particles) and in the low-density, 
high-temperature regime, the second virial coefficient gives the leading contribution to the 
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deviation of the equation of state from the non-interacting_ result of rewriting the 



59 



60]. 



grand canonical partition function as a cluster expansion 

The virial coefficient ([3]) turns out to be a simple, periodic (with period 2) but non- 
analytic function of the statistical parameter a, showing cusps in correspondence of all its 
bosonic points a — 2j. This quantity has been evaluated by different methods: one of them 
consists in an hard-disk-type regularization of the two-anyonic spectrum while another one 
is a based on path-integral approach yielding the two-body partition function, carried on 
by identifying the Lagrangian of the system with the one relative to the Bohm-Aharonov 
effect 6jJ. Eq. ((3]) is also retrieved by heat kernel methods, i.e. discretizing the two-particle 
spectrum through the introduction of a harmonic regulator potential and then directly con- 
sidering the problem in the continuum 62[. Finally, another method to get Eq. (3) is to 
use a semi-classical method, which nevertheless produces the exact quantum result 3a, 0|- 
Exact results for higher virial coefficients are not known, but a fair amount of information 
is av ailab ,e both f 0r the third v ir ial coefficient and fo r higher virial coefficients Q. 

The expression ([3]) is the exact quantum result for the hard-core case, corresponding to 
impose the vanishing of the two-anyonic wavefunctions in the coincident points (the limit 
configurations for which the coordinates of two anyons coincide). However, any arbitrary 
boundary condition for the wave-function is in principle admissible: it is the comparison with 
results from the microscopic interacting Hamiltonian that should fix the relevant boundary 
conditions to be imposed. The second virial coefficient for Abelian anyons in this general 



case has been studied in 



47, m 



63]. 



By relaxing the regularity requirements, allowing wave-functions to diverge for vanishing 
relative distance r between the anyons according to the method of self-adjoint extensions, 
it is possible to obtain a one-parameter family of boundary conditions. The hard-core limit 
corresponds to scale-invariance in a field theoretical approach 



46 



64 



67| , where the scale can 



be precisely related to the parameter of hard-coreness that will be defined below. The study 
of B2 shows that the results for hard-core case are rather peculiar: for instance, the cusps 
at the bosonic points of the hard-core case are a special feature of the scale-invariant limit, 
which is however absent for all the soft-core cases. On the contrary, cusps are generated 
at all the fermionic points for all the wavefunction boundary conditions, except just for the 
hard-core case. 

The relative two-body Hamiltonian for a free system of anyons with statistical parameter 



a, written in the bosonic description, is of the form 



H 



rel 



(5) 



where A = (A 1 , A 2 ) and A* 



1,2 and e u is the completely antisymmetric 



tensor). The corresponding single-particle partition function of the relative dynamics is 
Z re i = Tre~ l3Hrel , where /3 = 1/ksT. In order to proceed with the choice of a given self- 
adjoint extension, one has to define the space over which the trace above is performed. If we 
consider the radial component R n of the relative wave-function ip, the Schrodinger equation 
takes the form 



1 

M 



Id d (n + a) 
r dr dr r 2 



R n (r) = ER n (r) = ^R n (r) 



(6) 



with n even (choosing the bosonic description). Without any loss of generality, the statistical 
parameter can be chosen as a G [—1, 1]. Eq. (jBJ) is the Bessel equation and its general 
solution is given in terms of the Bessel functions 



R n {r) = AJ\ (kr) + BJ_\ n+a \(kr) . 



(7) 



For n 7^ the constant B must vanish in order to satisfy the normalization of the relative 
wave-function, while for n = (s-wave) arbitrary constants A, B are allowed. This yields a 
one-parameter family of boundary conditions for the s-wave solution: 



Rn(r) 



(const.] 



J\ a \(kr) 



' k 
I - 

K 



2\a\ 



J-\ a \{kr) 



(8) 



where a = ±1 and k is a scale introduced by the boundary condition. 
We will refer to 

e 



13k 2 



M 



(9) 



as the hard-coreness parameter of the gas. For e — > oo with a = +1 we retrieve the hard-core 
case (/0(O) = 0). If a — — 1, in addition to the solution (JS]), there is a bound state with 
energy Eb = —eksT = —k 2 /M and wavefunction 



Rn(r) 



(const.) K\ a \(nr) . 



(10) 



By proceeding as in 36(, and observing that only the s-wave energy spectrum is modified 
with respect to the hard-core case, one gets that the second virial coefficient for a generic 
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soft-core is given by 

B S 2 C \T) = BkHT) - 2\ 2 T I e-P EB 9(-a) + lim V 

R— >oo — ' 



s=0 



-RZs. 



e p m 



where 6{x) is the Heaviside step function, ko tS R is the s-th zero of J\ a \(kR) = 0, k s R is the 
s-th zero of (JSLj and B\ c - is the hard -core result ([3]). It is possible to rewrite Eq. (Till in an 
48(] as 



integral form 



B s 2 c -{T) = B h 2 c -{T)-2X 2 T { 
For e ^> 1 one gets 



Off 
7T 



Sin7TQ! 



o 1 + 2crcos7ra + t 2 l Q l 



(12) 



B S 2 C \T) = B^ C -(T) - 2X 2 T [e e 6(-a" 1 ' 



7T 



-I" 



sin 7r a + 



while for e <C 1 



^• C -(T) = B 2 C -(T) - 2X 2 T \a\ (1 - (re |a| + ■ • 



(13) 



(14) 



Near the bosonic point a = one has for £ ^ 0: 

1 



i?r-(T) 



4 



+ 2v(e)6(-a) 



\l + 0(a 2 



where i/(e) is the Neumann function defined by 

"°° cite* 



f e 



/o r(t + i) • 

From this expression one sees that B 2 C '(T) is a smooth function of the statistical parameter 
near a = 0. On the contrary, near the fermionic point |a| = 1 one has 



B S 2 C -(T) 



;(1-I«l) 2 + 



2e 



X 2 T + f a (£)(l-\a\) + 



where (er = ±1) are functions of e [not reported here], so that in general the soft-core 

case presents a cusp at \a\ = 1. 

The virial coefficient _Bf' c ' is plotted in FigJT]for some values of the hard-coreness param- 
eter £. The plot of the second virial coefficient clearly exhibits its smoothing in the bosonic 
points and its sharpening in the fermionic ones, as soon as the hard-core condition is relaxed. 
We observe that the restriction of B 2 C ' over the interval [0, 1] is not a monotonic function 
of a for each e. This not-monotonicity happens around £ ~ 1.4 in a narrow range of values 
of £ (1.344 < £ < 1.526), as emphasized in Figj2j 
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a l 



FIG. 1: B2 (in units of vs. the statistical parameter a for different values of the hard-coreness 
parameter e for Abelian anyons. The 6 upper curves are obtained for a = 1: from top to the 
bottom e takes the values 00 (hard-core), 10, 1.4, 1, 0.1 and 0. The curve below all those is 
obtained for a = — 1 and e = 0.1 and it has a shifted value at the bosonic points (the dotted lines 
just denote the x and y axes). We remind that B2 is periodic in a with period 2, so we plot it in 
the interval [—1,1]. B2 is symmetric with respect to all the integer value of a: we see from the 
figure the difference in the position of the cusps of the patterns, between the hard-core and the 
soft-core cases. 

B. Non- Abelian Anyons 



The main part of the present paper deals with the study of the low-density statistical 
mechanics properties of a two-dimensional gas of SU(2) non- Abelian Chern-Simons (NACS) 
spinless particles. The NACS particles are pointlike sources mutually interacting via a topo- 
logical non- Abelian Aharonov-Bohm effect 68[. These particles carry non- Abelian charges 
and non- Abelian magnetic fluxes, so that they acquire fractional spins and obey braid statis- 
tics as non- Abelian anyons. 

In order to proceed with the computation of the second virial coefficient of a free gas of 
NACS particles, we first introduce the NACS quantum mechanics 49_|, |69H72| considering 



50 



51] . The Hamiltonian describing the 



the general frame of soft-core NACS particles 
dynamics of the A^-body system of free NACS particles can be derived by a Lagrangian with 
a Chern-Simons term and a matter field coupled with the Chern-Simons gauge term j^ : 
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the resulting Hamiltonian reads 

- E w ( v *« v *» + V ^ v - Q ) ( 15 ) 



N 

H N 



0=1 

where M a is the mass of the a-th particles, V 2q = g|- and 

9 



dz a 2ttk a p z a - zp 

In Eq. (Tl5l) a = 1, . . . , N labels the particles, (x a , y a ) = (z a + z a , —i(z a — z a ))/2 are their 
spatial coordinates, and Q a, s are the isovector operators in a representation of isospin I. 
The quantum number I labels the irreducible representations of the group of the rotations 
induced by the coupling of the NACS particle matter field with the non-Abelian gauge field: 
as a consequence, the values of I are of course quantized and vary over all the integer and 
the half- integer numbers, with 1 = 1/2 being the smaller possible non-trivial value (/ = 
corresponds to a system of free bosons). As usual, a basis of isospin eigenstates can be labeled 
by I and the magnetic quantum number m (varying in the range — /, — / + 1, ■ ■ • ,1 — 1,1). 

The virial coefficients then depend in general on the value of the isospin quantum number 
I and on the coupling k (and of course on the temperature T). The quantity k in f fl5|) is a 
parameter of the theory. In order to enforce the gauge covariance of the theory the condition 
Awn = integer has to be satisfied. In the following we denote for simplicity the integer Ann 
by k: 

4ttk = k . (16) 

The physical meaning of k in the NACS model can be understood removing the interaction 
terms in by a similarity transformation: 

- 2 

H N — ► UHvU- 1 = Hfr = - E W 9 "" 9 " 

a a 

q H _^ uy H = y A (17) 



where U(z\, . . . , zn) satisfies the Knizhnik-Zamolodchikov (KZ) equation 73] 



k ~ w* g«^^J u{Zl ' ■ ■ ■ ' ZN) = ' (18) 

and tyij(zi, . . . , Zn) stands for the wavefunction of the iV-body system of the NACS particles 
in the holomorphic gauge. A comparison between the last equation and the KZ equation 
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satisfied by the Green's function in the conformal field theory shows that (4ttk — 2) cor- 
responds to the level of the underlying SU(2) current algebra. In 71j it is shown how 
^a(zi, . . . , zn) obeys the braid statistics due to the transformation function U(z±, . . . , zn), 
while ^h(zi, ■ ■ ■ , zn) satisfies ordinary statistics: ^a{zi, ■ ■ ■ , z^) can be then referred to as 
the NACS particle wavefunction in the anyon gauge. 

The statistical mechanics of the NACS particles can be studied by introducing the grand 
partition function S, defined in terms of the iV— body Hamiltonian and the fugacity v 
as 

oo 

E = ^v N Tre-f )H » . (19) 



7V=0 



In the low-density regime, a cluster expansion can be applied to H 59|, l60| : 



Z = exp^Vj2b n v n j , (20) 

where V is the volume of the gas (of course, for a two-dimensional gas V equals the area A) 
and b n is the n-th cluster integral, with 

l -Zu b 2 = Uz 2 -?f) (21) 



x A * A \ 2 t 
and = Tr e~^ HN being the iV-particle partition function. 

The virial expansion (i.e. the pressure expressed in powers of the density p = -j) is given 

as 

P = pk B T[l + B 2 (T)p + B 3 (T)p 2 + ...] , (22) 
where B n (T) is the n-th virial coefficient. The second virial coefficient B 2 {T) is written as 

B^T) = -% = A( l --^\ . (23) 



b\ \2 Z\j 

We assume that the NACS particles belong to the same isospin multiplet {\l,m >} with 
m = —I, The quantity Z\ = Tr e~^ Hx is then given by 

Z x = (21 + 1)A/X 2 T . (24) 



The computation of Z 2 = Tre~^ H2 is discussed in 55|, where the results for the hard-core 
case are presented. It is convenient to separate the center-of-mass and relative coordinates: 
defining Z = (z\ + z 2 )/2 and z = Z\ — z 2 one can write 

H 2 = H cm + H Tel = -—d z d 2 - -(V Z V- Z + V $ V Z ) , (25) 
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-0.25 - 



0.5 a 



FIG. 2: The narrow range of values of the hard-coreness parameter e within which the second virial 
coefficient B 2 (for Abelian anyons) is a non-monotonous function of a in the interval [0, 1]. The 
parameters are a = +1 for all the curves and e = 1.50, 1.45, 1.40, 1.35, 1.30 from top to bottom 
(-B2 is expressed in units of A T ). 

where fi = M/2 is the two-body reduced mass, V 2 = d 2 and 

V 2 = 9 2 + - . 



Q is a block-diagonal matrix given by 



2/ 



3=0 

with tjj = + 1) — 21(1 + 1)]. Z 2 can be then written as 



Z 2 — 2,A\j? Z 2 , 



(26) 



where Z' 2 = Tr rc ie @ H ™ X . The similarity transformation G(z,z) = exp { — ^ ln^z)}, acting 



as 



H re \ — > H' vcl — G 1 H rc \G, 
ty(z,z) — ► V'(z, z) = G~ x ^(z, z) 



(27) 



gives rise to an Hamiltonian H' rel manifestly Hermitian and leaves invariant Z 2 . The explicit 
expression for H' vA is 

^ = ~(v;v;+v;v;) , (28) 
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where V' z = 8 Z + fl/2z and V' 5 = d- z - Q/2z. 

By rewriting H' TcX in polar coordinates and projecting it onto the subspace of total isospin 
j, its correspondence with the Hamiltonian for (Abelian) anyons in the Coulomb gauge, 
having statistical parameter given by a s = Uj, becomes evident: 

2" 



8 2 18 1(8 

dr 2 r dr r 2 \ 89 



tU)j 



(29) 



The same analysis discussed in Section Hi Al shows that the radial factor of the j,j z — compo- 
nent of the relative (21 + l) 2 — vector wavefunction ip = e me R n (r) obeys the Bessel equation 



1 

M 



r dr dr 
whose general solution is 



Id d in + ujjY 
r— + ^ J - L 



R^(r)^ER^(r) = -R^(r) 



(30) 



R^(r) = A^J ln+UJj] (kr) + J-\ n+Wjl (kr) . (31) 

As already discussed in the previous Section III A\ B^ Jz can be nonzero only in the case n = 
(s— wave). Then the s— wave gives rise to a one-parameter family of boundary conditions 

k 



R, 



(const. 



J\^\(kr) 



a 



K 



2KI 



(32) 



where a = ±1, and Kj j x is a momentum scale introduced by the boundary condition. 
We refer to the {21 + l) 2 quantities 



M 



(33) 



as hard-coreness parameters of the system. The hard-core limit corresponds to Ejj z — > oo 
for all j,j z . 

We conclude this Section by observing that, according to the regularization used in 
621 ]. the second virial coefficient is defined as 



36 



B 2 (K,l,T)-Bt\l,T) 



2 A 2 ^ 



Z> 2 {K,l,T)-Z't L \hT) 



(34) 



(21 + l) 2 

where B^^'^i^T) is the virial coefficient for the system with particle isospin / and without 
interaction (k — > oo), which will be expressed in terms of the virial coefficients B^(T), 
B^(T) of the free Bose and Fermi systems with the considered general wavefunction bound- 
ary conditions. Furthermore, Z' 2 (k,1,T) — Z'^ 1 (l, T) is the (convergent) variation of the 
divergent partition function for the two-body relative Hamiltonian, between the interacting 
case in exam and the non- interacting limit (k — > oo). 
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III. SECOND VIRIAL COEFFICIENT 



In this Section we present our results for the second virial coefficient of a free gas of 
NACS particles: we will first study the hard-core case in Section IfflAl comparing in detail 
our findings with results available in literature 54j-l57|. We then study B 2 for non-Abelian 
anyons when general soft-core wavefunction boundary conditions, focusing the attention in 
particular to the isotropic boundary conditions. 



A. Hard-Core Case 



The hard-core case is obtained in the limit e 3 j z — » 00 for all j, j z . The second virial 
coefficient has been discussed in literature, and different results for B 2 have been presented 
56]: the differences between such results have been discussed, see in particular Ref. j^] 



54456 



and the comment 



56]. Our findings differ from results presented in 54H56| : in this Section, 



as well as in Appendices lATlBl a detailed comparison with such available results will be 
presented. 

For hard-core boundary conditions of the relative two-anyonic vectorial wavefunction, the 
quantity B 2 '^ entering Eq. (|34|) is found to be 55] 



1 21 



(21 



j=0 



1 + (-1V+ 21 



Bi(T) 



.1)3+21 



Bi(T) 



(35) 



Using for the hard-core case the known values B^(T) = —B 2 (T) = —\\? 



35j, one then 



obtains 



B^ L \l,T) 



1 



4 21 + 1 



(36) 



35 



To proceed further, we introduce a regularizing harmonic potential V 
whose effect is to make discrete the spectrum of H'y Using the notations of [2] (see pg.48), 
the spectrum consists of the following two classes: = e(2n + 1 + 7^) with degeneracy 
(n + 1), and E^ 1 = e(2n + 1 — 7^) with degeneracy n, where n is a non-negative integer and 
7j = ujj mod 2. It follows that the regularized partition function reads 



Z' 2 (K,l,T)-Z' 2 inA -\l,T) 



£(2j + 1) lin^ +l o ) {Z&j) - Z'M} 4 
3=0 



£->0 
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+ 



'-1)3+21 



{Z' e [( 7j + l)mod2]-Z' e (l)} 



(37) 



with 



n=0 



M2n+1+7 .) peQn+l-yM = 1 COsh [fo( 7 j - 1)] 

J 2 sinh 2 /3e 



The final result for the NACS gas in the hard-core limit is then the following: 



By c 



21 



2(21 + 



3=0 



1)3+21 



A 



T 



(7 -2 7 ,) + 



4 2Z 
1 - 



-1)3+% 



[( 7j - + l) mod 2 -If 



(31 



2 2 
Eq. (1551) is the main result of this Subsection: the dependence of B 2 on fc for some fixed 
Vs, and vice versa on / for some fixed fc's is represented in FigsJHB one sees from FigJH] a 
non-monotonic behavior of B 2 as a function of A;. In FigJUone can see that the values of B 2 
vs. I form two different groups, depending on integer and half- integer values of /. 

As a first check of Eq. (I38|) . we observe that the value of By 1 ' 1 ^ (l,T) of the free case 
(corresponding to the limit — > 0) is correctly reproduced: indeed, for a given I one 

0+ 



has in this limit Uj -+ O 1 * 1 , jj -+ 



Ti 



2 7i ->■ 0", [( 7i + 1) mod 2 - l] 2 -+ 0+, as it 



might. A more detailed discussion on the limit l/4irn — > is presented in Appendix IE1 

To further compare with available results, we observe that in 55| and 57J it was stated 
that the factors (— 1) 2/ should not appear in the expression of the two-particle partition 
function, or in the expression of the second virial coefficient (/ denoting the isospin quantum 
number of each particle). However, as pointed out in 56J and as further motivated in the 
following, such factors are needed. To clarify this issue it is convenient to make reference 
to the properties of the Clebsch-Gordan coefficients. The two (spinless) particles in exam 
have both isospin I and total isospin j: the Clebsch-Gordan coefficients express the change 
of basis, in the two-body isospin space, between the basis labeled by the individual magnetic 
isospin numbers mi, m 2 and the basis labeled by the total and magnetic isospin j, my. 

i i 

\lljmy) = (lmilm 2 \jmy)\lmilm 2 ) , 

mi=—l m,2=—l 

where the Clebsch-Gordan coefficients fulfill the symmetry property 



(lmilm 2 \jmj) = (— I) 21 J (lm 2 lmi\jmy) 
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100 



FIG. 3: F>2 vs. k = 4ttk for an hard-core NACS gas with 1 = 1/2 (blue open circles), I = 1 (red 
stars), 1 = 3/2 (green squares). In this and the following figures is in units of A 2 - (furthermore 
the line is just a guide for the eye, since k assumes only integer values). The values of B2 in the 
limit k — y 00 are given by — |, — y^, — yg and are correctly reproduced. 



Notice that the real spin of the particles is not taken in account in this consideration. With 
respect to the exchange of all the quantum numbers, the isospin two-body wavefunction 
corresponding to a state of total isospin j takes a factor (— l) 2/_jf = (— 1)-? +2 ' (being j an 
integer), so that the factor projects over the states for which the partition function 

can be evaluated in the bosonic basis, the factor — projects over those for which the 

partition function can be evaluated in the fermionic basis, from which Eqs. fl35|) . ( 136]) and 
(158]) can be obtained. 



Our result fl38l) differs also from the results presented in 54J , where a method of computa- 
tion of the second virial coefficient for hard-core NACS_gas based on the idea of averaging over 



all the isospin states is proposed. In particular, in 



54| the special cases / = 1/2, / = 1, and 



the large-K limit for two particles belonging to a representation / with lim^oo 



a < 1 



were considered. In the last limit the sum over all the resulting total isospins r < 21 is 
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approximated by an integral. The results are given by Eqs. (35), (36) and (38) of 541 ]: 



B%- c -{k,l = 1/2, T) = X 2 T 



B h 2 c - (k,l = 1,T) = X 2 T 



1 

1 



3 

4k 



3 
8F 



lim B% c - 

l— KX3 



A 2 ^ 



(M,T) 

while our corresponding results are 



B h 2 c - (k,l = l/2,T) = X 2 T 



20 

9k 



3k 2 



VA; > 1 



(39) 



- + a 



a 
~3 



3 

8P 



^• c -(A;,Z = l,r) = A 



T 



lim B\ c - (k, /, T) 



A 



T 



1 

12 

a 



3 

^8k 
14 8 
9A? ~ 3P 



VA; > 2 
VA; > 4 , 



(40) 



The limits I — > ao in (|39|) and ( 140]) are taken together with lim^ 00 / 2 /A; = a, where a is 
kept fixed and a < 1 in (139]) and a < 1/2 in f|4*0|) . The derivation from Eq. f|38|) of the 
three special cases above is reported in the Appendix |A] We notice that the asymptotic 
value found for B2 in the third case is expected to vanish for a = 0, while on the contrary 



this does not occur for the results of 



54] : indeed, a = corresponds to consider the limit 



Bf A \l,T) 



-1 A 



2Z+1 4 



which vanishes in the large /-limit. 



The difference between the results of Ref. 



54j and ours stands in a different averaging: 



while in 54j the virial coefficients are expressed as averages of the virial coefficients over the 
(2/ + l) 2 two-body states of isospin, in our case we take into account the effect of the isospin 
symmetry factor (— 1) J+2Z characterizing the states of total isospin j. 

Notice that Eq. (138]) for the second virial coefficient can be recovered using the approach 
presented in j^, as shown in Appendix |B] Indeed, one can find from Eqs. fl34]) - fl36|)-f[38l) fsee 
Appendix [B] for details) that the following expression for B%' c ~ given in Eq. (2) of 56| holds: 

1 21 



(21 



3=0 



where .B-f ,F (u;, T) is given by 



B, 



1 -1 + 45 - 25 2 , N even (odd) 



25 2 



N odd (even) 



(41) 



(42) 
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FIG. 4: i?2 vs. Z for an hard-core NACS gas with k = 1 (blue open circles), k = 2 (red stars), fc = 3 
(green squares). I varies over all the integer and the half-integer numbers: in the upper (lower) 
part of the figure the plotted values of B2 correspond to integer (half-integer) values of I. 



(u = N + 5 and N an integer such that < 8 < 1). The computation presented in Appendix 
IB1 shows that Eq. (2) of 56] is a correct starting point to study B\ c -\ however, notice that 
Eq. (3) of [56[ should be replaced with Eq. ( 1B3I) given in Appendix El 

We finally discuss in more detail the non-interacting limit I/Attk — > in order to clarify 
the meaning of Eq. (|36|) . In the limit k — > 00 the covariant derivatives in (fl5|) trivialize 
and the isospin becomes just a symmetry of the Hamiltonian, resulting in a pure (isospin) 
degeneration g = 21 + 1. Let H be the grand partition function, e the generic single- 
particle energy level for an assigned spectral discretization, and using the upper/lower signs 
respectively for the g-degenerate bosonic/fermionic single-particle states: notice that any 
value of g is allowed both in the bosonic and the fermionic case, since the statistics is not 
constrained by isospin. The expressions for the grand partition function, the pressure and 
the density are 



E(z,A,T) = H(lT 



ze 



MT9 



PA 



lnH(2,A,T)=Tj> ln (lT 



ze 



-fa 



and 



N = z-^\nE(z,A,T)=gJ2 



z-ie/fc =f 1 
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It follows 

( 



^ = ^ to(1 ^U* = * B TX;(£)V (43) 

9 n —r\ \9 J 



where B° denotes the n-th virial coefficient for spinless boson (/spinless fermion) without 
either isospin degeneration. Hence, denoting by B n the n-th virial coefficient in presence of 
isospin freedom one has 



OO oo 

P = k B T {P^tY—Bu = k B T P nB r, 

n=0 9 n=0 



and therefore 



therefore B 2 = 21+1^21 that is exactly what is written in Eq. (I36p . The result (136]) can be 
also understood by observing that all the virial coefficients for a system of NACS defined 
over a representation of isospin I tend, in the non- interacting limit k — > 00, to (— l) 2 ' times 
those of an ideal gas of identical /-spin ordinary quantum particles. In particular, the second 
virial coefficient of an ideal system of quantum s-spin particles is indeed: 

*^ = + TW- (44) 

in agreement with ( 1361) . The issue becomes much more complex for flux-carrying particles 
(finite k) having a non-zero spin, as discussed in 74-l76l| : however, for the true ideal spinor 



case a = and s = 1/2 it is B 2 (s = 1/2, T) = |75|, again in agreement with ( jH]) and 
the related ([36]). 

We conclude this Section by observing that a semiclassical computation of the second 
virial coefficient for a system of hard-core NACS particles reproduces Eq. ( 1381) : we remind 
that for an Abelian hard-core gas the semiclassical approximation {2), [33] yields the exact 
quantum result of {sgJ for B 2 . By extending such a computation to the hard-core NACS gas 
we find exactly Eq. (|38|) (details are not reported here). We mention that in literature it 
has been conjectured that the semiclassical approximation could give the exact expressions 
for all the virial coefficients in presence of hard-core boundary conditions 35] : the rationale 
for this conjecture is that for hard-core boundary conditions there are no other length scale 
besides At- Therefore, having established the extension to the non-Abelian hard-core case 
of the semiclassical computation of B 2 , one could in the future obtain information about 
higher virial coefficients for the hard-core NACS gas. However, we alert the reader that the 



21 



presence of other relevant length scales (other than At) in general prevent the semiclassical 



approximation from being exact: an explicit example is given in [35] . We conclude that 
for the soft-core NACS (that we are going to treat in the next Section) the semiclassical 
approximation is not expected to give the correct results. 



B. General Soft-Core Case 



If one removes the hard-core boundary condition for the relative (2/ + l) 2 -component two- 
anyonic wavefunction and fixes an arbitrary external potential as a spectral regularizator, 
then the spectrum of each projected Hamiltonian operator ifj can be represented as the 
union of the spectra of (2j + 1) scalar Schrodinger operators, one for each j 2 -component, 
endowed with its respective hard-coreness parameter Ejj z (as shown in Appendix [C]). As 
discussed in Section [II Bl one then ends up with a set of (21 + l) 2 (in principle independent) 
parameters £jj z , which are needed to fix the boundary behavior. They can be organized in 
a (21 + 1) x (21 + 1) matrix: 



/ 



£ 0,0 


£1,1 


£2,2 


£1,-1 


£1,0 


£2,1 


£2,-2 


£2,-1 


£2,0 



£21+1,21 



(45) 



\ £22+1,-22-1 £21+1,-21 £22+1,0 / 

Proceeding as in the previous Subsection IIII At one has then for the general soft-core NACS 
gas the following expression for the second virial coefficient: 

21 i r-i , I 1\ 14-22 1 I _]\j+2l 



(2l 



j=0 jx=-j 



1 _|_ (_l)j+2« I 
O B 2 (UjjTjEjj,) + — 



where _B.f (u;^, T, e^J is the soft-core expression entering Eq. ( [121) : 



r poo 

Bffa^Ejj,) = ^' c -(5 j7 T)-2A 2 \ e £ ^6(-a) + ^- (sinvr^ 



-B 2 (u)j,T, £j,j z 
(46) 



1 + 2a(cos7r5 j ) +t 2 l^l 

(47) 
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with Sj = (ujj + l)mod2 — 1, and B^iojj, T, Sj j z ) is the previous expression evaluated for 



ujj — > Uj + 1: 



Bl {ah T, tu .) = i**(r„r)-2Aj (e-e(-<x) + ^ (sin^ , + (|r ,, + (2|rjl 

(48) 

with Tj = Ujmod2 — 1. Eq. ( )46|) is the desired result for a NACS ideal gas with general 
soft-core boundary conditions. 

To perform explicit computations, we consider in the following the simple case in which 
the isotropy of the hard-coreness parameter is assumed within each shell with assigned 
isospin quantum number I. In other words, £jj z = Ej and the matrix fj45 jl then reads 

/ 



-J,Jz 



Sq Ei 
E X Ei 



£21+1 \ 
£21+1 



(49) 



\ £21+1 £21+1 ■ ■ • £21+1 / 

When all the element of the matrix (09]) are equal, we will use the notation £jj z = e. In 
such a completely isotropic case, Eq. ( |46|) takes the simplified form 

2/ 



B s 2 c -(k,1,T) 



1 



(2/+l) : 



^(2j + l)5f(^,T,e) 



j=0 



where 



1 + (-\y +21 



mod 2 



(50) 



(51) 



In Figs. [Mini we show, for three values of the isotropic hard-coreness parameter e, the depen- 
dence of i?2' c ' on k for some fixed Z's, and vice versa on I for some fixed fc's. FigJT] evidences 
that for suitable values of e the second virial coefficient may change sign and have strong 
variations. From Eq. fl50|) it is possible to see that the values of -B|' c ' (ft, /, T) corresponding 
to semi-integer / and k = 1 are independent of I, depending only on e and T (see Figs J6|8lfT0|) . 
In fact Eq. (I5T|) yields that for I semi-integer (/ = 1/2, 3/2, ■ 
and therefore 



and k = 1 one has Vj = ±| 



k = 1, 1 = n H — , T 
2 



(52) 



(with n = 0,1,2, 
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FIG. 5: -B|' c '(/c, I, T) vs. k for a soft-core NACS gas with / = 1/2 (blue open circles), I = 1 (red 
stars) and 1 = 3/2 (green squares): in all e = 0.1 (k varies over the positive integers). 



For I = 1/2, i.e. the lowest possible value of I for non-Abelian anyons, the assumption of 
isotropy (e ,o — £ o an -d £i, m = e i with m = 1, 0, —1) yields: 



B S 2 C - 



k, I = -, T 
2' 



-B B (k>i,T, £1) 



As example, let us consider the case I = l/2,4nn = 3: 



k 



3,/ = -,T 



o 



"2' T ' £ ° 



(53) 



(54) 



(similar results can be found for other values of k). 

For the isotropic soft-core system, special boundary conditions are the s-channel of soft- 
coreness (for which the p-wave is assumed to be hard-core, e± = 00), and the p-channel of 
soft-coreness (for which, vice versa, the s-wave is assumed to be hard-core, Eq = 00). In 
order to assure the physical soundness of the virial expansion, ksT has to be much higher 



47l | than the energy of the eventual bound state Eb associated to the wavefunction (11 01) . 



Hence, for both these channels the virial expansion is meaningful provided that we take 
a = +1 in Eq. f )12p . and the virial coefficients for these two channels are 



B S 2 C - 



and 



b s 2 c - 



k = 3,1 



3,7 = 5 



n as, r. 24 r° . e -"H- 1 '' 



s— channel 



T 

24 



7T 



dt- 



p— channel 



24 I vr 



dt 



g— eit^— 5/6 

1 + v^^ + t 1 / 3 



(55) 



(56) 
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FIG. 6: B 2 c '{k, I, T) vs. I for k = 1 (blue open circles), k = 2 (red stars) and fc = 3 (green squares), 
with e = 0.1 (/ varies over the integer and the half-integer numbers). 

The previous equation clearly shows that the depletion of B 2 with respect to the hard-core 
value — 24 is the result of the anyonic collisions allowed by the soft-core conditions. If 
the four parameters of the whole matrix are taken to be identical Eq = E\ = e ("complete 
isotropy" of the hard-coreness matrix), the expression for the virial coefficient reduces to 

B t> r*= 3 ,/= d + i r* e -« /6 r " 2 



2' y 24 i ttJo V 1 + t 1 + v^t^ + t 1 / 3 

(57) 

The dependence of this quantity on e becomes more evident by representing the e variable in 
logarithmic scale, as shown in FigJTTl The hard-core limit value B^'^/X^ = —1/24 predicted 
by (|38p is asymptotically approached, although for extremely high e: e.g. for e = 10 17 it is 
B\- c ' j\\ w —0.05, which deviates from the asymptotic value by a ~ 20%. We conclude that 
even an extremely small soft-coreness may have a significant impact on B 2 and therefore 
on the thermodynamical properties. At variance, for small values of e, the extension of 
the analysis presented in 48j for soft-core Abelian anyons allows to compute the value of 
B 2 for the limit case e = 0: for e = one gets B s 2 - C - (k = 3, 1 = |,T) j\\ = -13/24. The 
monotonically increasing behaviour of B 2 in e is evident from (I57I) . and consistent with an 
approach towards an hard-core (hence more repulsive) condition. 



25 




20 k 30 



FIG. 7: B2 C '(k,l,T) vs. for I = 1/2 (blue open circles), I = 1 (red stars) and I = 3/2 (green 
squares), with s = 
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FIG. 8: -B|' c '(/c, /, T) vs. Z for k = 1 (blue open circles), k = 2 (red stars) and fc = 3 (green squares), 
with e = 1.4. 

IV. OTHER THERMODYNAMIC AL PROPERTIES 



In this Section we remind the virial expansions for the main thermodynamical quantities 
of the system studied in this paper, namely the gas of NACS particles endowed with general 
boundary conditions for the wavefunctions and we discuss how these quantities read at the 
order pA^ of the virial coefficient, in order to highlight the role played by the second virial 
coefficient B 2 computed in the previous Section. 
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FIG. 9: -B|' c '(fc, I, T) vs. k for I = 1/2 (blue open circles), I = 1 (red stars) and I = 3/2 (green 
squares), with e = 10. 

The thermodynamical quantities are associated to the virial coefficients {B n (T)} of the 
equation of state (featuring the expansion of the pressure in powers of the number density 
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77j one has the following virial 



p). As discussed in statistical mechanics textbooks 
expansions for the pressure P, the Helmholtz free energy Ah, the Gibbs free energy G, the 
entropy S, the internal energy E and the enthalpy H (A being the area) : 

PA 



Pressure : 

Helmholtz free energy : 
Gibbs free energy : 



Nk B T 

A H 
Nk B T 

G 



k>l 



lo 



g (px 2 T ) - 1 + \ B k+i 



P k ; 



fe>i 



Nk B T 



k+i p ; 



k>l 



S 1 d 

— = 2 - log(pA^) - ^ - _ (T B k+1 ) p k 



Nk B 
Internal energy : 
H 



k>l 



E 



Nk B T 



H ( 

Enthalpy : = 2 + ( Bk+1 

B k>i ^ 



fc dT 

k>l 

_l T dB k+1 



P 



Nk B T ' ^ \ ^ k dT 

Using the previous expressions, stopping at the lowest order of the virial coefficient p 
(i.e. pAy) and using the fact that for a general NACS ideal gas one has 52 (T) oc T" 1 , one 
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FIG. 10: B 2 c '(k,l,T) vs. £ for = 1 (blue open circles), k = 2 (red stars) and = 3 (green 
squares), with s = 10. 

can obtain the thermodynamical quantities at the lowest order of the virial expansion: in 
particular we find 

PA 

1 + B 2 p ; 



Nk B T 



A 



H 



Nk B T 
G 
Nk B T 



E 



Nk B T 
H 



log(pA^) - 1 + B 2 p 
--\og{p\ 2 T ) + 2B 2 p ■ 

1 + B 2 p ; 

2 + 2B 2 p 



Nk B T 

(at the lowest order of virial expansion, the entropy and the heat capacity at constant volume 
do not depend on B 2 ). Using the expression of B 2 given by Eq. ( )46l) one can obtain the 
deviation of the various thermodynamical quantities from their ideal gas value. 

An important consequence of the previous results is that at the order pA^ one finds 



E = PA : 



(58) 



Eq. ( 1581) is an exact identity for 2D Bose and Fermi ideal gases, valid at all the orders of 



the virial expansion Similarly, for the soft-core NACS ideal gas, at the order pA^, one 



has H = 2E, which is also exact at all orders for 2D Bose and Fermi ideal gases [2|. Further 
investigations on the higher virial coefficients are needed to ascertain if the equation of state 
fl58|) is exact (at all orders) for a general soft-core NACS ideal gas. 
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FIG. 11: -B|' c ' as a function of the hard-coreness parameter e in logarithmic scale, for k = 3 and 
1 = 1/2, in the completely isotropic case (eo,o = £ i,m = s, for m = 1, 0, —1). 

V. CONCLUSIONS 

In this paper we studied the thermodynamical properties of an ideal gas of non-Abelian 
Chern-Simons particles considering the effect of general soft-core boundary conditions for 
the two-body wavefunction at zero distance: we determined and studied the second virial 
coefficient as a function of the coupling k and the (iso)spin I for generic hard-coreness 
parameters. A discussion of the comparison of obtained findings with available results 
in literature for systems of non-Abelian hard-core Chern-Simons particles has been also 
supplied. We also found that a semiclassical computation of the second virial coefficient 
for hard-core non-Abelian Chern-Simons particles gives the correct result, extending in this 
way the corresponding result for Abelian hard-core anyons. We have also wrote down the 
expressions for the thermodynamical quantities at the lower order of the virial expansion, 
finding that at this order the relation between the internal energy and the pressure is the 
same found (exactly) for 2D Bose and Fermi ideal gases. Further studies on the higher virial 
coefficients are needed to establish the eventual validity of the obtained relation between 
the pressure and the internal energy for a general soft-core NACS ideal gas. 
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Appendix A: Computation of virial coefficients in special cases 



To perform the comparison with Ref. 54| . we compute in the following the virial 
coefficients for a NACS gas in the hard-core limit in the special cases considered in 54J: 
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It is w = — |, Wi = — |, w 2 = I and the cu/s are such that \ujj\ < 1: therefore, by usinj 
l2"l). one has 
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of even 2/ (the opposite one is similar). We define j cr u as the maximum integer such that 
Uj < 0, and x cr u = \/2l: it can be verified that \tUj\ < 1 for all j. We have 
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Appendix B: Comparison with Ref. [56] 



With the notation used in the main text, Eq. (2) of Ref. 56J for NACS particles in the 
hard-core limit reads 
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with u = N + S and N an integer such that < 8 < 1. Eq. ( 1B1I) can be derived as in the 
following: with the notation jj = Ujmod2, using Eqs. (I34j) - (136j) - (l38j) one has 
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that is nothing else than the (IBlj) itself. 

Notice that Eq. (3) of [56] should be replaced with 
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where the upper and lower signs refer to the cases of even and odd iVj's: Eq. (IB3|) is 
equivalent to our formula ( 138|) and to Eq. (IBlj) . 

We discuss now the limit a = — — )■ : we observe that by a careful inspection it is 
possible to conclude that Eq. (30) of js5| and Eq. (3) of 5f| do not tend in this limit to 
the correct value B^ 1 (l,T) given in Eq. ( )35i) . However, the manipulation of the corrected 
version (lB3j) above presented reproduces (as expected) the value B\ c - — > — ^^zTT ^ or tt ~ ^ 0- 
Indeed for vanishing coupling constant a, using the same convention used above (upper and 
lower choices referring to the cases of even and odd Nj respectively), one has 
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as it might. 
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Appendix C: Spectrum in the soft-core non-Abelian case 



In Subsection IIIIBI we stated that the spectrum of the mult i- component projected Hamil- 
tonian operator H'j can be represented in general as the union of the (2j + 1) spectra of the 
corresponding scalar Schrodinger operators. That follows from the following remark: the 
non-Abelian generalization of the soft-core expression (jHJ) is Eq. ( 132]) . By denoting 
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the non-Abelian generalization of ([6]) is the (21 



-dimensional matricial equation 



( A ( r ) 



V ■■• 



o 

A(r) 







A(r) 



\ 



( i? °'° \ 



/ 



R 



1,-1 



R 



1.0 



V 



/ R° ° \ 



E 



J 



R, 



1,-1 



pl,0 



\ 



( R° '° \ 



M 



J 



R 



1,-1 



R 



1.0 



V 



(CI) 



/ 



The hard-disk regularization i?Q-' z (r = R) = for all j,j z discretizes the energy spectrum; 
let then Spjj,(R) be the (discretized) spectrum of the component equation A(r)R 3 Q Z = 
ERq- ,z restricted over the domain in which Rq 3z has hard-coreness parameter £j,j z , and 
Sp(R) be the (discretized) spectrum of Eq. (ICip in the domain in which any component 
R Q has the respective assigned hard-coreness parameter £j>j'. If Ejijr G Spji^ z (R) for 
some (j',j' z ), then also Ej/jt G Sp(R), because 
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(C2) 



and all the null components of the vector trivially fulfill whichever hard-coreness conditions, 
in particular the assigned sequence {£jj,} G {[0, oo)}^ 2 '" 1 " 1 -' 2 . In conclusion, the spectrum in 
the non-Abelian case can be written as the above-mentioned union of spectra, which will 
automatically include all the possible relevant energy degenerations to be considered in the 
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partition function for the computation of the virial coefficients. 
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